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Spin-entanglement of two electrons occupying two spatial regions - domains - is expressed in a compact form in 
terms of spin-spin correlation functions. The power of the formalism is demonstrated on several examples ranging 
from generation of entanglement by scattering of two electrons to the entanglement of a pair of qubits represented 
by a double quantum dot coupled to leads. In the latter case the collapse of entanglement due to the Kondo effect is 
analyzed. 
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1 Introduction Based on the peculiar behaviour of en- 
tangled quantum states, in 1935 Einstein, Podolsky and 
Rosen argued that quantum mechanical description of real- 
ity is not complete yj. Today this article is Einstein's most 
cited publication, but quantum entanglement is not consid- 
ered a paradox. In fact, the ability to establish entangle- 
ment between quantum particles in a controlled manner is 
a crucial ingredient of any quantum information process- 
ing system 121. Also, the study of entanglement provides 
insight into the nature of many-body states in the vicinity 
of crossovers between various regimes or points of quan- 
tum phase transition fS) . 

In realistic hardware designed for quantum information 
processing, several criteria for qubits (DiVincenzo's check- 
list) must be fulfilled ||4]: the existence of multiple identifi- 
able qubits, the ability to initialize and manipulate qubits, 
small decoherence, and the ability to measure qubits, i.e., 
to determine the outcome of computation. Quite generally, 
the parts of an interacting system are to some extent entan- 
gled. However, fully entangled qubit pairs are required for 
such applications [jj]. This leads to the question of how to 
quantify the entanglement. 

The entanglement of binary quantum objects (qubits) 
can be quantified by the entanglement of formation, a no- 
tion which for pure states reduces to the von Neumann en- 



tropy ||5]|6]|7l. This measure is convenient since it is pro- 
portional to the number of fully entangled pairs needed to 
produce a given entangled state. Also, it is monotonically 
related to the concurrence, which can be evaluated by the 
Wootters formula JD as a simple function of the density 
matrix. 

A possible realization of a qubit pair can be achieved by 
confining electrons to two spatial regions. In this case one 
has the freedom to choose between spin or charge of the 
electron to represent a qubit. The entanglement measures 
must account for both possibilities. Additional complica- 
tions arise due to the indistinguishability of the particles 
and possible states of multiple occupancy ll9l [T0lfm[T2l . 

On the experimental side, it seems that among several 
proposals for quantum information processing systems, the 
criteria for scalable qubits can be met in soUd state struc- 
tures consisting of coupled quantum dots II13II14I . The abil- 
ity to precisely control the number of electrons ifTsl and 
the evidence for spin entangled states have been reported 
in GaAs based heterostructures II16II17I . It has also been 
demonstrated that in double quantum dot systems coher- 
ent qubit manipulation and projective readout are possible 

EH. 

Here we explicitely evaluate the concurrence in terms 
of spin-spin correlation functions between the two subsys- 
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cupying a separate quantum dot A and B, can be described 
in the basis of single spin-i states s =| or J, as \^ab) = 
X^ss' oiss'\s)a\s')b- Such a state, Fig.lita), is factorizable 
I^'ab) = \(j)A)\4>B) if and only if a^ia^^ = a^a^ as 
easily checked by explicit construction. For general coeffi- 
cients ass' the state is not factorizable and its entanglement 
can be quantified by concurrence C as introduced by Hill 
and Wootters ||6l, 



C = 2\a 



Ti"iT 



(1) 



Consider now a more general problem of two electrons 
in a state where the system cannot directly be reduced to an 
equivalent system with (pseudo) spin degrees of freedom 
only. Let us take for example two electrons on a lattice |[T9l 
described by a state 



(d) 
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Figure 1 (Color online) (a) Two electrons (qubits) local- 
ized in separate quantum dots A and B. (b) Two delocal- 
ized electrons in separate measurement domains A and B; 
the probability of finding one electron in a domain is equal 
for A and B: ua = ns = 1. (c) Scattering of two electrons, 
(d) Double quantum dot system coupled to the leads: the 
chemical potential is set such that each of the dots is occu- 
pied by one interacting electron on average. 



tems. In this manner we easily account for entanglement 
between subsystems containing delocalized electrons. We 
demonstrate the convenience of our approach by evaluat- 
ing the time-dependence of concurrence in a toy example 
of two electrons in a three site chain and in a more realis- 
tic entangler based on scattering of electrons in a quantum 
wire. 

On the other hand, entangled qubits can also be ex- 
tracted from interacting many-body systems. Our formu- 
las are readily generalized also for this case and even for 
finite temperature. As an example, we discuss the case of 
two quantum dots coupled serially to metallic electrodes. 
In this case, entanglement of a pair of electrons that are 
confined in a double quantum dot may collapse at low tem- 
peratures due to the Kondo effect even for a very weak 
coupling to the leads. 

The paper is organized as follows. Sec. |2] introduces 
the entanglement measure for two delocalized electrons, 
which can be simplified for some special cases elaborated 
in Sec. [3] In Sec. |4]we illustrate the convenience of our 
approach with numerical examples mentioned above. Con- 
clusions are given in Sec. |5] 

2 Entanglement of two delocalized electrons A 

pure state of two spins represented by electrons, each oc- 



\^) 



ss' 



IS js' 



|0), 



(2) 



where cj^ creates an electron with spin s on site i on the lat- 
tice with the total number of sites N . Sites of the lattice are 
assumed to be ordered as for example in Fig. [lib) and (i, j) 
corresponds to the summation over all sites i ~ 1, ...N and 

« < J ESI- 

Such states may arise when two initially unentangled 
electrons in wave packets approach each other, then inter- 
act and finally become again well separated in distinct re- 
gions A and B, Fig. [Tic), where they could be extracted 
for further purposes. Alternatively, such states can be real- 
ized in various correlated electron systems and, in order to 
study them theoretically, elaborate many-body techniques 
are sometimes needed. Moreover, usually not the state it- 
self but only the correlation functions are available. 

Therefore it is advantageous to express the entangle- 
ment in terms of spin-spin correlation functions. The spin 
operators for a single electron occupying domain A (or B) 
are expressed as the sum of operators for sites i within the 
given domain. 



'^A{B) 



E 

ieA{B) ss' 



(3) 



where (T^^, are the Pauli matrices. 

The concurrence is given by the eigenvalues of the non- 
Hermitian matrix pp, where p and p are the reduced density 
matrix and its time reverse, respectively JS]. For axially 
symmetric problems, where both (S'lS'^j.^^jS') = and 

{^31^^)1"^) = as well as mS^AiB)S%^A)\1') = and 
(^\'^A(B)^B{A) l'^) ^ be expressed in a compact 

form 011 



C = max(0,CTx,q|), 
=21(5+5^)1-2 



q, =2|(5+5+)|-2(pX), 



(4) 
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where S+^g^ = {S^^g^Y ^ T.^eA(B) are spin rais- 

ing operators for domains A (or B) and 

Pmb) = 'i^Ul-n,^-s) (5) 

are spin-s projectors operating in domains A (or B) with 
nis = cJ^Cis. Fermionic expectation values required in 
Eq. (|4|i are then given as 

(5i5+)=Ev5*<, (6) 

M 

M 

where [ij] in Eqs. ^ corresponds to the summation over 
all pairs i,j such that i G A and j G B. Concurrence 
formula Eq. (HI is valid as long as double occupancy of 
sites is negligible, {riisUis') — > 0. It is assumed that the 
wave function is normalized, 

(^^1*')= J2 \^tf^ = i- (7) 

We stress that the correlation functions in Eq. ^ can be 
evaluated for pure or mixed states, the latter arising, for 
instance, due to the finite temperature or by tracing out the 
environment degrees of freedom, e.g., leads in Fig. 1(d). 

3 Special cases In states with the SU(2) symmetry 
{S^S%) = {S\Sl) = {S%SI^) and the concurrence for- 
mula Eq. dUi simplifies further to a function depending on 
only one spin invariant (Sa • Sb), 

Cab = max (^0, -2(Sa • Sb) - . (8) 

The concurrence is expected to be significant whenever en- 
hanced spin-spin correlations indicate A-B singlet forma- 
tion. 

If \^) is an eigenstate of the total spin projection S^^^ 
the concurrence is given solely with the overlap between 
I If-) and the AB-spin-flipped state |#) = \S^Sg>F). If 
5jot = ±1 the concurrence is zero, while for S^^^ = 

C = CTi=2|^^1]>Jj|, (9) 

M 

which is a generalization of the concurrence formula Eq. ([T]i 
to N sites. 

The concurrence formulas, Eq. (HI, remain essentially 
the same if the state corresponds to the system in con- 
tinuum space, i ^ r = {x,y, z), the only change being 
integrations over the corresponding measurement domains, 

C=\f /(<^o-<^i)*(</^o+<^i)d'rid3r2|, (10) 
JaJb 



where (ps = </'s(ri,r2) = (ri, r2; S'j'^') are singlet and 
triplet amplitudes for 5 = and 5=1, respectively. 

Another interesting special case is the wave function 
which is a linear combination of entangled Bell AB- 
pairs, 

4 

13=1 [ij] [tj] 

where for each pair of sites one can introduce the Bell 
basis \ij,(3) 13, 



U, 1) = 






*J,2) = 






■ij, 3) = 


— (c^ 




U,4) = 







(13) 

In this case, the concurrence is given with a simple expres- 
sion C = | 62|_ 

4 Numerical examples Here we use concurrence for- 
mulas in practice. We evaluate the concurrence for a few 
examples of interacting electrons on a lattice described by 
the following generic hamiltonian 

H = - '^iUjcl^Cjs +h.c.) + E Uijn^srijs'. (14) 

ijs ijss' 

For simplicity we take the electron-electron interaction con- 
stant up to some distance, i.e., Uij = J2m=o 

4.1 Two qubits on three sites First we consider two 
electrons on three sites, where site i = 1 corresponds to 
the measuring domain A and sites i = 2, 3 to the domain 
B. We take the Hubbard model, M = 0, and two non-zero 
hopping matrix elements, ti2 and t23. In the limit i23 > tri 
a bonding orbital is formed between the sites 2 and 3 and in 
the ground state of the system there is a single electron in 
each of the domains. For large U, the ground state is a spin 
singlet formed between site i = 1 and the bonding orbital, 
with the excitation energy J to the triplet state. 

If in each of the domains there is precisely one electron 
and the state is an eigenstate of total spin projection, S^^^ = 
0, Eq. (01 simplifies to 

C = 2|ViV<2+WiVf3l- (15) 

Let us put the electrons to the system in an initially sep- 
arable state consisting of a spin up electron in A and the 
other electron with spin down in the bonding orbital of B, 
\]]/{0)) = cL4=(4i + C31 Because the inital state is 
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Figure 2 (Color online) Full line represents concurrence 
C{t) for a three sites system with ^12/^23 = 1/10' 
U/t23 = 5 and Jeff/i23 = 0.023 (singlet-triplet energy 
difference). Dashed line corresponds to the strong coupling 
limit results Ceff (^) with the same Jeff. Domains are doubly 
occupied with the probabihty P2. 



composed of different energy eigen-states, the Rabi oscil- 
lations occur. In the strong coupling limit, U, ^23 ^ tri the 
system is described by the Heisenberg model with antifer- 
romagnetic coupling Jeff ^ ^ii2l{U — 2i23) between the 
site 1 and the bonding orbital. In this limit the Rabi oscilla- 
tions occur due to the singlet-triplet splitting and the time 
evolution of concurrence is given by Ceff = | sin Jsfft\. 

For generic values of parameters additional states, for 
example the states when the site 1 is doubly occupied, be- 
come relevant. We show the concurrence for such a case 
in Fig.|2]and compare it to the simplified expression given 
above. The simple behaviour is partly reproduced but ad- 
ditional oscilations arise on other time scales. We plot also 
the probability that both electrons simultaneously occupy 
one of the domains, P2 = 1— (ni|+ni^)+2(ni|ni|). The 
oscillations of this quantity occur on the time scale given 
by the characteristic time of tunneling events, 

4.2 Two flying qubits Next we consider two flying 
qubits, i.e., two electrons on an infinite one-dimensional 
lattice with the hamiltonian Eq. ( fT4] i and tij = io for j = 
Tobe specific, let one electron with spin f be confined 
initially to region A(i ^ ~L) and the other electron with 
opposite spin to region B {i ^ L), Fig. [3] The simplest 
initial state consists of two wave packets with vanishing 
momentum uncertainties Ak 0, with momenta fc > 
and q < for the left and the right wave packet, respec- 
tively. After the collision, the electrons move apart with a 
probability amplitude Ai^^{k,q) for non-spin-flip scatter- 
ing and a spin-flip amplitude A^^ (fc, q). Concurrence after 
the collision is then readily expressed from Eq. ^ as 

C = 2\Mn{k,q)Mu{k,q)\. (16) 




(=180/(„ C=0.97 



-120 -80 -40 40 80 120 

sites i 




Figure 3 (Color online) (a) Time evolution of spin f and 
I electron density during the scattering for the Hubbard 
model (M — 0) with U = to and Ak — 7r/50 and fco near 
the concurrence maximum, (b) C as a function of fco: (i) 
the Hubbard model (M = 0) with U = to and Ak = 0; 
(ii) M = 3: with Ak = (full line), and (iii) Ak = tt/IO 
(dashed). 

Note that C = 1 when non-spin-flip and spin-flip ampli- 
tudes coincide in accord with recent analysis of flying and 
static qubits entanglement Il21ll22l|23ll24l . 

More general initial wave packets with finite Ak are 
defined with appropriate momentum amplitudes (f>^{k) and 
4>l{q) for spin f and |, respectively. The concurrence as 
follows from Eq. (|9]l 

C = 2\ 11 MJ^{k,q)Mu{k,q)\h{k)<Pi{q)\'dkdq\ 

(17) 

consists of a coherent superposition of scattering ampli- 
tudes 111 9]. 

The simplest example is the Hubbard model where the 
scattering amplitudes can be obtained analytically for the 
case of one-dimension, A^||(fc, g) = 1 + A^||(fc, g) = 
(sin fc - sin g)/ [sin fc - sin g + lU/ (2io)] EH. In Fig.Ha) 
the time evolution of spin t and j electron densities is 
presented. In Fig. |3jb) the corresponding concurrence for 
wave packets with a well defined momentum fcg for U = to 
is shown, together with a longer range interaction case, 
M = 3, for a sharp momentum (full line) and for a Gaus- 
sian initial amplitude 0x(g) — (t)^{~k) with Ak — tt/IO 
(dashed line). An interesting observation here is a substan- 
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Figure 4 (Color online) (a) Time evolution of spin t and 
I electron density for U = 2to, M = 5 with fco = O.Stt, 
Z\fc| = Ak^ = 7r/20. (b) C{t) for Gaussian packets with 
various fco, M = 5, U ^ 2to and Ak = 7r/20. At t = 
the separation between the packets is 2L = 10/ Ak. 



tial reduction of the concurrence due to coherent averaging 
in Eq. ivh . Additionally, electrons will be completely en- 
tangled at some kinetic energy comparable with the repul- 
sion, U ^ 2to(l — cos fco), where non-spin-flip and spin- 
flip amplitudes coincide. In Fig. HJa) and Fig. |5ja) some 
representative additional examples of time evolution of in- 
teracting wave packets are shown. 

The concurrence formula Eq. (|9]l is derived for elec- 
tronic states when double occupancy is negligible, which 
in our case is strictly fulfilled only asymptotically when 
the electrons are far apart. However, Eq. (|9]l can be eval- 
uated at any time t and the resulting C{t) can serve as a 
measure of entanglement during the transition from initial 
to final state. In Fig.Utb) and Fig.|5jb) C{t) corresponding 
to parameters of Fig.|4ja) and Fig.|5ja) is shown. Concur- 
rence oscillations can be interpreted as a response to the 
finite time duration of electron-electron interaction - ex- 
change - where the model can be approximately mapped 
onto an effective Heisenberg model as in the case of three 
sites presented above. 

4.3 Qubit pairs in coupled quantum dots One of 

the simplest realizations of a solid state qubit is a serially 
coupled double quantum (DQD). We model such a DQD 



Figure 5 (Color online) As in Fig.|4] with parameters U = 
2to, Af = 5 but different fco = OAn, Ak^ = tt/IO and 
Aki - tt/5. 

using the two-impurity Anderson Hamiltonian 

H= {erii + Un^^nii) + VnAnB (18) 

i=A,B 

s 

where cl^ creates an electron with spin s in the dot i = A 
or i = B and Uis = cl,,Cis is the number operator The 
on-site energies e and the Hubbard repulsion U are taken 
equal for both dots. The dots are coupled to the left and 
right noninteracting tight-binding leads with the chemical 
potential set to the middle of the band of width 4to- Each of 
the dots is coupled to the adjacent lead by hopping t2 and 
the corresponding hybridization width is = (^2)^/^0- 
Schematically this setup is presented in Fig.llJd)- The dots 
are additionally coupled capacitively by a inter-dot repul- 
sion term VuAn-B- 

When a DQD is attached to leads the low tempera- 
ture physics is to a large extent the same as that of the 
two-impurity Kondo problem |i26j. There two impurities 
form either two Kondo singlets with delocalized electrons 
or bind into a local spin-singlet state which is virtually 
decoupled from delocalized electrons. The crossover be- 
tween the regimes is determined by the relative values of 
the exchange magnetic energy J and twice the Kondo con- 
densation energy, of order the Kondo temperature given 
by the Haldane formula, Tk = ^?7i"/2 cxp(— 7re(e + 
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Figure 6 (a) Spin-spin correlation for V/U = 
0,1/3,2/3,1 (full lines) and V/U = 5/4 (dotted) for 
U/r = 12, r/ta = 0.1. (b) Probabilities for parallel 
(lower curves - dashed) and anti-parallel (upper curves - 
full) spins of electrons in the DQD for V/U ratios as in 
(a). Note that the probability for parallel spins for V/U ^ 
5/4 is almost zero (dashed-dotted), while P|| < 1/2 for 
J/Tk < 1000 (dotted); the probabiUties do not sum to 1. 
The deficiency (which goes to zero as U ^ oo) is due to 
the states with double occupancy on at least one of the dots. 

U)/2r). The competition between extended Kondo and 
local-singlet states occurs rather generally in systems of 
coupled quantum impurities II27II28I[29]| . 

A qubit pair represented by two electrons in a DQD and 
in the contact with the leads acting as a fermionic bath can 
not be described by a pure state. In the case of mixed states 
of qubit pairs as studied here the concurrence is related to 
the reduced density matrix of the DQD subsystem as in the 
previous case of delocalized electrons pairs. If tn2)/U is 
not small the electrons can fluctuate between the dots (and 
to the leads) and such charge fluctuations introduce addi- 
tional states with zero or double occupancy of individual 
dots |.9,,U ,| . As pointed out by Zanardi in the case of simple 
Hubbard dimer lITTl the entanglement is not related only to 
spin but also to charge degrees of freedom which emerge 
when repulsion between electrons is weak or moderate. 

For systems with strong electron-electron repulsion con- 
sidered here, charge fluctuations are suppressed and the 
states with single occupancy - the spin-qubits - dominate 
and the concept of spin-entanglement quantified with con- 
currence can be applied. We again use spin-projected den- 
sity matrix and consider only entanglement corresponding 
to spin degrees of freedom. Due to doubly (or zero) oc- 
cupied states arising from charge fluctuation on the dots 
(caused by tunneling between the dots A and B or due to 




4 8 12 

UIT 



100 


r(b) 


<S;S,>=-l/4 ^ 






/ 




: [//r=i2 


/ 


10 




/A. 








1 





0.2 0.4 0.6 0.8 1 1.2 

VIU 



Figure 7 (a) Charge fluctuations of the domain A, An\ = 
{n\) - {uAf (short-dashed), (S^ • Sb) = -1/4 (long- 
dashed) and C = (full) in the [U/T, J/Tk) plane, (b) 
(Sa • Sb) = -1/4 (long-dashed) and C = (full) in the 
{V/U, J/Tk) plane. 

the exchange with the electrons in the leads), the reduced 
density matrix has to be renormalized. The probability that 
at the measurement of entanglement there is precisely one 
electron on each of the dots is less than unity, Pn < 1, and 
the spin-concurrence is then given with 

Cii = C/Pii, (19) 
where Pn = P^^ + i^,, and P^^ = {pIp^ + PjP^), 

parallel spin alignment, respectively. Such a procedure cor- 
responds to the measurement apparatus which would only 
discern spins and ignore all cases whenever no electron, or 
a electron pair would appear at one of the detectors at sites 
AorB. 

Here we present results for the zero temperature limit 
obtained by the projection operator method as in Refs. 1301 
I3ni32ll33ll34l . The temperature dependence of the concur- 
rence for the case of = is given in Ref . |j 3? !l . Expecta- 
tion values (...) in the concurrence formula Eq. ^ are now 
calculated using the ground state therefore {Sj^S^) — 
andC|| < 0. We consider the particle-hole symmetric point 
with n = 2 and e + U/2 + = 0. 

Qualitatively, the concurrence is significant whenever 
enhanced spin-spin correlations indicate inter-dot singlet 
formation. As shown in Fig.|6la) for U/P = 12 and P/tg = 
0.1, the correlation function (Sa • Sb) tends to —3/4 for J 
large enough to suppress the formation of Kondo singlets, 
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but still J /U ^ 1, that local charge fluctuations are suf- 
ficiently suppressed. In particular, the local dot-dot singlet 
is formed whenever singlet-triplet splitting superexchange 
energy J > Jc ^ 2Tk- With increasing V ^ U, and 
above U, the probability for singly occupied spin states, 
Pii = P^i + P|| is significantly reduced, Fig.|6jb), which 
also leads to reduced spin-spin correlation. 

The entanglement between qubits quantified by con- 
currence is small in the regime where the Kondo effect 
determines the ground state. The Kondo screening trans- 
fers the entanglement between localized electrons to the 
mutual entanglement of localized and the conducting elec- 
trons [[361 • For V ^ U the Kondo temperature is enhanced 
and the corresponding Kondo ground state is competitive 
towards the localized singlet state. In Fig. |2la) we present 
a phase diagram of regions with finite and zero entangle- 
ment, separated approximately by the J ^ 2Tk line. In 
the strong coupling limit the boundary coincides with the 
(Sa • Sb) = -1/4 dashed fine. For V > U + Tk the 
Kondo screening is inhibited and the concurrence is in- 
creased, Fig.|2tb). However, the probability for doubly oc- 
cupied sites is not small in this regime and the concept of 
spin-qubits is not appropriate there ifSTl . 

5 Conclusion In this work we showed that the appli- 
cability of the concurrence, an entanglement measure orig- 
inally restricted to two distinguishable particles (spins), can 
be extended to two regions - measurement domains - oc- 
cupied by two indistinguishable particles (electrons). The 
proposed approach allows for the analysis of entanglement 
in a variety of realistic problems, from scattering of fly- 
ing and static qubits, the former being represented as wave 
packets with a finite energy resolution, to the time evolu- 
tion of static qubits due to electron-electron interaction or 
externally applied fields. A generalization to systems de- 
scribed as a mixed state and containing more than two elec- 
trons is possible. As an example, a double quantum dot 
system occupied, on average, by two interacting electrons 
is presented. 

We acknowledge support from the Slovenian Research 
Agency under contract Pl-0044. 
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